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Abstract 

Stereoscopic arrays of Imaging Atmospheric Cherenkov Telescopes allow to reconstruct gamma- 
ray-induced showers in 3 dimensions, which offers several advantages: direct access to the 
shower parameters in space and straightforward calorimetric measurement of the incident energy. 
In addition, correlations between the different images of the same shower are taken into account. 
An analysis method based on a simple 3D-model of electromagnetic showers was recently im- 
plemented in the framework of the H.E.S.S. experiment. In the present article, the method is 
completed by an additional quality criterion, which reduces the background contamination by a 
factor of about 2 in the case of extended sources, while keeping gamma-ray efficiency at a high 
level. On the other hand, the dramatic flares of the blazar PKS 2155-304 in July 2006, which 
provided H.E.S.S. data with an almost pure gamma-ray sample, offered the unique opportunity 
of a precision test of the 3D-reconstruction method as well as of the H.E.S.S. simulations used 
in its calibration. An agreement at a few percent level is found between data and simulations for 
the distributions of all 3D shower parameters. 

Key words: Gamma-ray astronomy, H.E.S.S., Stereoscopy, Cherenkov telescopes, 
3D-reconstruction, Analysis method 
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1. Introduction 

During the last four years, an important gain in sensitivity was achieved in the very-high- 
energy domain of gamma-ray astronomy, resulting in an increase of the number of confirmed 
sources by a factor of about ten. This progress was made possible by the new generation of 
Imaging Atmospheric Cherenkov Telescopes, and particularly by stereoscopic arrays of such de- 
tectors. Presently and in the near future, the number of individual telescopes in the different 
arrays is limited: 5 in H.E.S.S.-II d], 4 in CANGAROO-3 ||2|], H.E.S.S.-I U and VERITAS 
101 and 2 in MAGIC-II f3\. In the long term, much larger arrays are considered in order to 
improve both the sensitivity (through the rejection of hadronic showers and the increase of the 
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effective detection area) and the angular resolution. The 3D-reconstruction of gamma-ray show- 
ers developed by M. Lemoine-Goumard et al. in the framework of the H.E.S.S. experiment is 
particularly well suited to large arrays since correlations between different stereoscopic views of 
the same shower are taken into account. Furthermore, in this method, gamma-rays are selected 
on the basis of physical properties: the rotational symmetry of an electromagnetic shower and 
the lateral spread of the corresponding Cherenkov photosphere at shower maximum. The total 
number of Cherenkov photons emitted by the shower is then reconstructed, yielding the gamma- 
ray energy in a straightforward way. This 3D-model was successfully applied in some H.E.S.S. 
studies, both on point-like sources (e.g. the blazar H2356-309 |9]) and on extended sources (e.g. 
the supernova remnant RX J0852. 0-4622 |lOfl, also named Vela Junior). 

In the present article, we first describe a significant improvement in the model which has 
been achieved since the original publication |6j]. In the 3D-reconstruction method, the shower 
parameters are obtained by a maximum likelihood fit to the charge contents of the pixels of each 
camera, with the single constraint of rotational symmetrjQ. However, in the version described 
in 10], no goodness-of-fit parameters were used in the further selection of gamma-ray showers. 
This conservative strategy, in which only a loose convergence criterion was required, was due to 
the fact that light fluctuations at the pixel level were not perfectly controlled and that correlations 
between the contents of different pixels were not taken into account. Moreover, in order to be 
reliable, a goodness-of-fit parameter based on individual pixel charges would require a complete 
modelling of the effects of the night-sky background. In section |2] below, we show that the sit- 
uation is simpler at the level of each triggered telescope, and that we can directly compare the 
total photo-electron yield per camera predicted by the model to the actually measured one. Fur- 
thermore, in non-triggered telescopes, the model should predict a modest light yield for genuine 
gamma-ray showers. Such comparisons allow to define self-consistency conditions which should 
be satisfied for most electromagnetic showers but not necessarily for hadronic ones. Applying 
this additional quality criterion to H.E.S.S. data on an extended sourceH, the hadronic background 
remaining after using standard 3D-cuts is reduced by a factor of about 2, while gamma-ray effi- 
ciency is kept at a high level (~ 80%). The significance of a weak gamma-ray signal can thus be 
improved by a factor 1.24 with respect to the standard 3D-analysis. 



The performance of the 3D-analysis in the H.E.S.S. experiment was estimated from gamma- 
ray simulations which were also used to calibrate the energy measurement in this method. The 
simulation procedure used in this study consists of two parts: (a) development of gamma-ray- 
induced showers in the atmosphere (program KASCADE 1 14]) and (b) complete description of 
the detection chain (program SMASH 11511 ). namely optics, triggering electronics and readout 
electronics of the H.E.S.S. telescopes. In step (b) the average evolution of the detector with time 
is taken into account, since the efficiency of the whole acquisition chain (mirror reflectivity, light 
collectors efficiency, quantum efficiency of the phototubes, electronic charge measurement) is 
regularly monitored on the basis of ring-like images produced by muons jl ill il2ll when tele- 
scopes are triggered individually. The dramatic flares of the blazar PKS 2155-304 in July 2006 
1 1311 in which the average intensity of the source was about 7 times that of the Crab nebula. 



'The effect of the weak geomagnetic field (20/jT on the H.E.S.S. site, due to the South Atlantic Anomaly) is taken 
into account in the simulations and in the calculation of the performance of the method. 

^In the case of an extended source, no cut on the angle between the source nominal position and the shower axis is 
appHed. 
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Figure 1: Left panel: distribution of the variable u = ln(Q3o/2real)/'"' gamma-ray showers simulated at a zenith angle 
o/ 18°. From top to bottom: showers viewed by 2 telescopes (black), 4 telescopes (green) and 3 telescopes (red). Right 
panel: distribution of the variable u = ln(Q3i) / Qxnl) for gamma-ray showers simulated at different zenith angles f and 
viewed by 2 telescopes. From top to bottom: ( = 18°, 37°, 48° and 53°. Gamma-rays were simulated according to a 
power-law energy spectrum with a spectral index of 2.2. 



provided us with a gamma-ray beam with a very small contamination of hadronic background, 
offering a unique opportunity of checking both the performance of the 3D-model and the quality 
of the whole simulation chain with a high accuracy. The corresponding tests are described in 
section |3] 



2. Characterising the quality of the likelihood fit 

The 3D reconstruction method is based on a likelihood fit in which the shower parameters 
in space are obtained from the charge contents (in singlephoto-electron units, hereafter p.e.) of 
individual camera pixels (see formula (1) of reference [6]). Images were previously submitted 
to a cleaning procedure in order to remove isolated clusters of pixels with a small charge which 
likely result from the night sky background. The likelihood fit was restricted to those pixels re- 
tained by the cleaning procedure and to their immediate cleared neighbours, as explained in fl. 
On the other hand, starting from the fitted shower parameters, it is possible to calculate the total 
photo-electron yield Qt,d expected from the 3D model in each camera. This quantity is also cal- 
culated for telescopes which were not triggered by the processed event and for which no image 
is available, since, in the framework of the H.E.S.S. trigger and acquisition system |7], they are 
not read out. It should be emphasised that, for a given triggered telescope, the predicted charge 
as well as the measured charge Qreai to which it will be compared, encompass the whole 
field of view (except for those pixels invalidated by the calibration) and not only the restricted 
sample of pixels used in the fit. Extending the comparison to the full field of view allows a better 
rejection of hadronic showers, generally more patchy than those of electromagnetic ones, and in 
which light is spilled over a region larger than the image core. 
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Figure 2: Distributions of the variable u = ln(23D/2ical)/o'' gamma-ray showers simulated at a zenith angle o/ 18° and 
viewed by 2, 3 and 4 telescopes (from left to right). Dotted lines show the distributions fitted according to formula Q. 



2.1. Quality criterion on images from triggered telescopes 

The compatibility between the model and the image can be characterised by the variable 
u - ln(23D/6reai)- Simulations of gamma-ray-induced showers at different energies and zenith 
angles were used to find the expected u distributions. In these simulations, the complete detection 
chain is taken into account. Simulated images were further submitted to the usual analysis proce- 
dure based on the 3D-model. The distribution of the variable u is shown in Fig.[T](left panel) for 
gamma-ray showers simulated at a zenith angle of 18°, according to a power-law energy spec- 
trum with a spectral index of 2.2. Using different simulated spectra (e.g. with spectral indices 
greater than 3) does not significantly change the distribution. This will be further illustrated in 
section[3]on real data from the blazar PKS 2155-304 whose spectral index is 3.3 (see Fig.[T3]). 
The distribution peaks at m = as expected and is slightly asymmetric as shown in Fig.[T] right 
panel, for events triggering 2 telescopes at different zenith angles. Fig.[T](left panel) also shows 
some dependence on the telescope multiplicity, i.e. the number nj of triggered telescopes. The 
difference between the case nj - 2 and that of higher multiplicities is related to the property of 
rotational symmetry of the shower which is assumed but not really checked if nj - 2, whereas 
it acts as an additional constraint for nj- > 3, resulting in a slightly broader distribution of the 
variable u in the latter case. 

For M < the distribution is well described by a Gaussian, whereas it behaves as an expo- 
nential one for large positive values of u. Therefore the probability density function piu) can be 
conveniently described by the convolution of a Gaussian distribution, with mean value uq and 
rm.s. s, and of an exponential one with mean value IjA, namely : 

p{u) ^ K e-^" fK(\{'^—!^\ with freq(x)=^= f e-^^'^du (1) 

the normalization factor K being given by : 

A^s^ 

K ^ A exp(Auo —) 

The parameters mq, s and A were fitted to the histograms obtained from simulations. The corre- 
sponding fitted distributions are indicated by the curves superimposed to the histograms in Fig.|2l 
The interest of the analytic parametrisation ([B is that the corresponding cumulative distribution 
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Figure 3: Distributions of the goodness variable P{u) for showers viewed by 2, 3 and 4 telescopes (from left to right). 
Top: gamma-ray showers simulated at a zenith angle o/ 18°, according to a power-law energy spectrum with a spectral 
index of 2.2. Bottom: hadronic showers as observed by H.E.S.S. at a zenith angle of 20° in a field of view with no 
significant gamma-ray emission. 



is easily calculated as : 

Q{u) = J" p{v) dv = freq ^ 

The distribution of Q{u) is thus expected to be uniform for gamma-ray showers, whereas hadronic 
ones will exhibit peaks for small values of Q{u) (corresponding to strongly negative values of u) 
and for large values (corresponding to strongly positive values of u). It is convenient to fold this 
distribution in order to obtain a variable whose behaviour is the same as a.x^ probability. This is 
done by using the following variable : 

P{u) = 2 min(e(M), 1 - Q{u)) (3) 

P{u) will be called a "goodness variable" in the following. The distributions of P{u) for those 
simulated gamma-ray showers referred to in Fig. [T] (left panel) are shown in Fig. [3] (top panels) 
and are fairly uniforrrH On the other hand, distributions of P{u) for hadronic showers detected 
by H.E.S.S. at a zenith angle of 20° in a field of view with no significant gamma-ray emission 
are shown in Fig.[3](bottom panels) for telescope multiplicities nj- = 2, 3 and 4 and show a strong 
accumulation at very low values. 

In practice, the dependence of the parameters mq, s and A as functions of the telescope mul- 
tiplicity tij and the zenith angle ^ is determined on the basis of simulations. The parameter mq 



A 



(2) 



'of course, these distributions need not be exactly uniform, the effect of the cuts on P(u) being determined by 
simulations. 
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Figure 4: Gamma-ray showers simulated at a zenith angle o/ 18°; the goodness variables P{u) of the same event viewed 
by two telescopes are plotted versus each other 



is always found to be compatible with within the statistical errors, which confirms the quality 
of the reconstruction method. As far as the other two parameters s and A are concerned, the 
empirical formulae shown in Table [TJ obtained from simulations, give a satisfactory description 
of the u distributions in all conditions of observation. 





s 


A 


2 


0.183 -0.186 cos^ + 0.134 cos^^ 


min(6.3, 9 


358 -4.411 cos 


3 


0.084 + 0.087 cos ^ 


7.141 


- 3.528 cos ^ 


4 


0.108 + 0.062 cos ^ 


8.703 


-4.986 cos ^ 



Table 1: Formulae giving the parameters s and A defined in equation Q as fimctions of the telescope multiplicity nj and 
of the zenith angle (. These formulae were obtained from simulations. 

For each event, formulae ([TJ, Q and (|3]l and those given in Table [T] allow to calculate one 
goodness variable P{u) per triggered telescope. For a given event, these quantities are correlated ; 
in particular, in the case of genuine gamma-ray showers, there is a clear positive correlation 
between the different values of P{u). Therefore, if P(u) is required to be greater than a minimal 
value pminfor all triggered telescopes, the gamma-ray selection efficiency of such a cut is greater 
than p"'^. . This is verified from Fig. |4]in which two values of P(u) of each gamma-ray shower 
are plotted versus each other. The regions close to the axes are depopulated with respect to that 
surrounding the diagonal. The slight excess close to the origin (about 20 events for a total of 
8400 simulated gamma-rays) is compatible with the contribution of hadron photo-production 
in the primary interaction (occurring with a probability of =s 2.8 x 10"^ ||171 ), which is taken 



Figure 5: Hadronic showers as observed by H.E.S.S. at a zenith angle of 20° in a field of view with no significant 
gamma-ray emission: the goodness variables P(u) of the same event viewed by two telescopes are plotted versus each 
other 



into account in the simulation. On the other hand, the same cut provides a powerful hadronic 
rejection, as shown in Fig. |5]in which two values of P{u) of each hadronic shower are plotted 
versus each other Those events for which at least one probability is very small accumulate in 
two "walls" and a spectacular spike is visible at the origin. 

2.2. Quality criterion on images from non-triggered telescopes 



nj 


In go 


go for ^ = 20° 


O-D 


2 


4.367-0.092 cos^ 


72 


0.109 + 0.405 cosf 


3 


max(4.5, 4.296 + 0.288 cos 


96 


0.184 + 0.259 cos f 



Table 2: Formulae giving Qq (in p.e.) and the standard deviation cro of the distributions o/ln^so as functions of the 
telescope multiplicity nj and on the z,enith angle These formulae were obtained from simulations. 

The preceding simulations of gamma-ray showers were also used to find the distribution of 
the light yield expected from the 3D-model in non-triggered telescopes. The same sample of 
simulated gamma-ray showers as in the preceding subsection was used to find the distributions 
of the charge Q^d predicted by the reconstruction method in non-triggered telescopes. They are 
shown in Fig. |6l both for nj - 2 (i.e. 2 non-triggered telescopes) and for nj - 3 (i.e. one 
non-triggered telescope). It can be verified that typical values of this undetected charge range 
between 50 and 100 p.e. Furthermore, it turns out that the distribution of \n{QT,r,) is practically 
Gaussian (see Fig.|2|. 



c 




Q3D [Pe] 

Figure 6: Distribution of the charge Q30 expected from the 3D-model in non-triggered telescopes for gamma-ray showers 
simulated at a zenith angle of 18°. The solid line corresponds to events with 2 non-triggered telescopes and the dotted 
line to events with only one non-triggered telescope. Gamma-rays were simulated according to a power-law energy 
spectrum with a spectral index of 2.2. 




ln(Q3D) ln(Q3D) 

Figure 7: Distributions of \n{Q->,o) for gamma-ray showers simulated at a zenith angle 0/ 18°. Left panel: events with 
2 non-triggered telescopes (nj = 2) and Q30 > 20 p.e. Right panel: events with 1 non-triggered telescope (nj = 3j. 
Gamma-rays were simulated according to a power-law energy spectrum with a spectral index of 2.2. 
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For each zenith angle and for a given value of «/•, the average value ln(2o) and the stan- 
dard deviation cr^j of this distribution are determined from simulations. This allows to define the 
normalised Gaussian variable v - ln(23£)/Go)/cD, whose square behaves as ax^ with one de- 
gree of freedom. It is thus possible to use the corresponding probability Piv^) as an additional 
goodness variable. There is one such variable per non-triggered telescope which, for gamma-ray 
showers, should be uniformly distributed between and 1. This is actually verified in Fig.[8](top 
panels), which shows these distributions for gamma-ray showers simulated at a zenith angle of 
18° for iij -2 and nj - 3. In the case of rij - 2, the small fraction 2%) of gamma-ray events 
whose predicted charge Qt,£, in a non-triggered telescope is less than 20 p.e. (see Fig.|6]l has been 
removed from the histogram on the top left panel of Fig. [8] These events, which correspond to 
showers falling at large distance from the centre of the array, have a very small probability Piy^) 
and would have been rejected by the standard cuts too. There are no such events for tij - 3. 
Conversely, Fig. [8] (bottom panels) shows the distributions of the probability Piy^) for hadronic 
showers as observed by H.E.S.S. at a zenith angle of 20° in a field of view with no significant 
gamma-ray emission. The accumulation at very low probabilities is very clear. The dependence 
of In 2o and that of the standard deviation ctq of the distributions of In Qo,^ on the telescope mul- 
tiplicity iij and on the zenith angle ^ are obtained from simulations and are given in Table|2] The 
values found for are comparable to typical cutoff values for the charge collected in a given 
telescope, as used in other H.E.S.S. analysis methods lfl6ll ifisll . 

2.3. Overall quality criterion and performance of the method 

For any event reconstructed by the 3D-model, 4 goodness variables are available, one for 
each telescope, whatever the trigger decision. For gamma-ray showers, these quantities should 
be uniformly distributed between and 1, whereas hadronic ones should most of the time lead to 
low values. The new quality criterion requires the goodness variables of all triggered telescopes 
to be higher than pnig - 0.05 and the goodness variables of all non-triggered telescopes to be 
higher than /?dark - 0.07 . This ensures that the data are compatible with the basic assumptions 
of the model (mainly rotational symmetry) and with the fitted values of the shower parameters. 
After this first step, the standard selection cuts described in are applied. They require the 
compatibility of the fitted parameters with those of an electromagnetic shower, namely : 

• compatibility between the depth of shower maximum and the number of Cherenkov pho- 
tons (formula (2) of reference |6]); 

• restriction on the shower 3D-width, in practice on the "reduced 3D-width", whose distri- 
bution is almost independent of the zenith angle (formula (3) of reference |@]). 

We have checked that the distributions of the different goodness variables defined above did not 
change significantly after applying the two preceding standard cuts. Finally, in the case of a 
point-like source, a cut is applied on the angle between the reconstructed direction and that of 
the source 10, namely 6 < 0.1° / cos 

The new method is compared to the standard 3D analysi^in Table[3]for an extended source 
(i.e. without any cut on the shower direction) and in Table H] for a point-like source. For each 
kind of analysis, gamma-ray selection efficiencies e are obtained from simulations assuming a 



Compared to the H.E.S.S. standard analysis, based on the mean-scaled widths and lengths of shower images, the 
standard 3D analysis provides a gain of 1.12 in significance LSI. 
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Figure 8: Distributions of the probability P(v^) obtained from non-triggered telescopes for events with 2 non-triggered 
telescopes (nr =2, left panels) and events with I non-triggered telescope (nj = 3, right panel). Top panels: gamma-ray 
showers simulated at a zenith angle of 18°, according to a power-law energy spectrum with a spectral index of 2.2. 
Bottom panels: hadronic showers as observed by H.E.S.S. at a zenith angle of 20° in a field of view with no significant 
gamma-ray emission. 



Zenith angle 


18° 


37° 


46° 


3D 




0.92 


0.91 


0.91 


standard 


Rh 


9.95 


9.69 


10.51 


analysis 


Gext 


2.89 


2.85 


2.94 


New 




0.78 


0.78 


0.77 


3D 


Rh 


21.01 


18.54 


18.05 


analysis 


2ext 


3.59 


3.34 


3.26 




Gainext 


1.24 


1.17 


1.11 



Table 3: New 3D analysis as compared to the standard one for an extended source (no angular cut is applied): eg is 
the gamma-ray selection efficiency, Ri, the hadronic rejection factor and Qext = % the quality factor. The gain in 
significance (Gainextj is given by the ratio of quality factors. 
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Zenith angle 


18° 


37° 


46° 


3D 




0.58 


0.60 


0.58 


standard 


Rh 


5292 


4233 


2966 


analysis 


Qpoint 


41.9 


38.71 


31.82 


New 




0.51 


0.52 


0.51 


3D 


Rh 


8468 


6124 


4029 


analysis 


(2point 


47.11 


40.92 


32.21 




Gain point 


1.12 


1.06 


1.01 



Table 4: New 3D analysis as compared to the standard one for a point-like source (angular cut 6 < 0. 1° / cosf j; Cg is 
the gamma-ray selection efficiency, Ri, the hadronic rejection factor and Qpoint = the quality factor The gain in 

significance (Gainpointj is given by the ratio of quality factors. 



3D standard 
analysis with 
«7- > 3 


Zenith angle 


18° 


33° 


51° 


Rh (extended) 
Rh (point-like) 


16.1 
6870 


13.9 
4980 


11.3 
2170 



Table 5: Hadronic rejection factors for extended and point-like sources in the 3D standard analysis in which the restric- 
tion nj- > 3 is applied. Results are reproduced from ^ which uses slightly different values of the zenith angle. For the 
sake of comparison with Tables\3\and^ linear interpolations in cos ( yield Rf, (extended) = 13.7 for ( = 37° and 12.1 
for f = 46°. Similarly. Rj, (point-like) = 4560 for ^ = 37° and 3050 for f = 46°. 



power-law energy spectrum with photon index 2.2. The loss of genuine gamma-ray events due 
to the new criteria (~ 15% on average) is almost energy-independent above 100 GeV. The corre- 
sponding rejection factors R for hadrons (i.e. the factors by which the numbers of hadronic events 
at the trigger level are reduced after applying selection cuts) are obtained from real H.E.S.S. data 
taken in fields of view with no significant gamma-ray emission. In the case of a weak source, the 
significance of the signal is proportional to the quality factor Q - e y/R which therefore charac- 
terises the sensitivity of the method. The ratio of the quality factors in the new 3D analysis and 
in the standard one respectively indicates the gain in significance obtained from the new method 
for a source detection. When values of p^ng and /^dark lower than those given above are used, 
the quality factors and the gain are found to be rather stable. The values chosen for these two 
probabilities have the advantage to significantly reduce the residual contamination of the selected 
sample. For an extended source at low zenith angles, the background is reduced by a factor of 
2, which is important for morphological as well as for spectral studies. The gain in sensitivity, 
namely 1.24 at low zenith angles for an extended source, corresponds to that of an observation 
time increased by 50% in the framework of the standard 3D analysis. For point-like sources at 
low zenith angles, the gain is lower albeit still significant and the background is reduced by a 
factor 1.6. It should be noted that the improvement of the method is stronger for events trigger- 
ing two telescopes which are the more numerous. The standard 3D analysis had sometimes been 
applied with the restriction nj- > 3 in order to achieve a higher sensitivity (see e.g. |0]). With the 
new analysis, the hadronic rejection factors obtained without any cut on the telescope multiplicity 
are greater than or comparable to those obtained in the standard 3D method with nj > 3. From 
Table |5] in which the corresponding factors are reproduced from reference |6], it can be verified 
that the rejection factors obtained in the new analysis with nj > 2 are higher by 20 to 30%. 
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3. Testing the method and the simulations with the almost pure gamma-ray sample from 
the flaring period of the blazar PKS 2155-304 in July 2006 

The exceptional VHE gamma-ray flare of PKS 2155-304 during the nights MJD 53944 and 
MJD 53946, in which the average intensity of the outburst was about 7 times the flux observed 
from the Crab Nebula 1 13], provided us with an unprecedented data-set taken under uniform ex- 
perimental conditions. The very high quality data featuring an abundance of gamma-rays with a 
very small contamination (less than 3%) of hadronic background effectively provides us with a 
gamma-ray test-beam . To correct for these remaining hadronic events, the reflected background 
model subtraction method |8] has been applied. In this case, the OFF region from which the 
background is estimated is the exact symmetric counterpart of the ON region with respect to the 
telescope axis, so that acceptances in both regions are the same. 



The simulations used for the comparison were custom-made to fit the actual experimental 
conditions during the observation period, such as the night sky background rate (40 MHz) and 
the optical efliciency of the telescopes, which is directly measured via images of muon rings 
II12II . The latter is important for the calibration of the energy measurement, which relies on the 
reconstruction of the number of Cherenkov photons in the shower and is therefore particularly 
sensitive to variations in the relative optical efficiency. The data-set has been divided into sam- 
ples corresponding to given intervals of zenith angle used in the simulation and the comparison 
between data and simulations was carried out for all these samples. The gamma-ray energy spec- 
trum was simulated according to a power law, using the spectral index of 3.3 fitted from the data 
of the two flaring nights. 



Observable quantities that are characteristic of the 3D-Model and related to crucial exper- 
imental parameters have been compared systematically, thereby examining the degree of ac- 
curacy of the match. In the following, all the comparisons refer to a zenith angle of 18° but 
they are representative for the whole range of zenith angles covered by this study. For each 
observable X, whose mean values are Xd (real data) and Xs (simulations), the relative bias 
2(Xd - Xs)/iXd + Xs) = AX/X is given in Table |6] All these quantities are at most of a few 
percent. The half-widths half-maximum (HWHMd for data, HWHM, for simulations) of all 
histograms shown in Fig. |9][l0][TT]and[T2](left panel) have been calculated. Table|6]shows that, 
for all observables, the differences A{HWHM) = HWHMd - HWHM,, are smaller than the bin 
width. 



3.1. Trigger-related quantities 

Fig. |9] shows the agreement for trigger-related quantities. The distribution of the number 
of pixels hit by photons in the camera is in very good agreement with the simulation. Espe- 
cially the rising edge of the distribution, which is particularly sensitive to the actual experimental 
conditions, such as the night sky background rate, has been reproduced by the simulation to a 
high degree of accuracy. A second trigger-sensitive variable is the amount of charge in the third 
most luminous pixel in each telescope. As the trigger condition requires a threshold of 4 photo- 
electrons in at least three pixels in each sector of the camera, the comparison of the charge in the 
third most luminous pixel allows us to scrutinise the accuracy of the trigger. 
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Figure 9: Comparison of observable quantities related to the trigger between the data (red points with error bars) and 
3D-Model simulations (black histogram). Left panel: Number of pixels per telescope. Right panel: Charge in the third 
most brightest pixel of one of the telescopes. 




Figure 10: Comparison of obsen'able quantities related to the atmospheric profile between the data (red points with 
error bars) and 3D-Model simulations (black histogram). Left panel: Altitude of the shower maximum. Right panel: 
Slant depth (grammage) of the shower maximum. 
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3.2. Atmospheric profile 

Other environmental conditions studied here include the atmospheric profile (desert atmo- 
sphere) used in the simulation. Whereas the altitude of shower maximum, as reconstructed by 
the 3D model, is a purely geometric parameter, the slant depth of shower maximum expressed 
in g/cm^ (grammage) is determined on the basis of the atmospheric profile. Since the distribu- 
tions of the first quantity show a good agreement between data and simulations (see Fig. [lOl 
left panel), any discrepancy between the true atmospheric profile and that of simulations would 
induce a difference between the distributions of the slant depth obtained from simulations and 
data respectively. This is however not observed in Fig. [TO] (right panel), which validates the 
atmospheric profile implemented in simulations. 

3.3. Shower characteristics 




5 10 15 20 25 30 35 40 45 50 0.0005 o.ool 0.001 s 0.002 0.0025 0.003 0.0035 0.004 0.004s 0.005 



3D-width reduced 3D-width 

Figure 1 1 : Comparison of observable quantities related to the atmospheric profile between the data ( red points with 
error bars) and 3D-Model simulations (hlaclc histogram). Left panel: 3D-width. Right panel: Reduced 3D-width. 



Even more pivotal than the agreement of environmental parameters is the accordance of the 
observable quantities upon which gamma-rays are separated from hadronic showers. In the 3D- 
model, this separation relies on the reduced 3D-width of the observed shower. Noting that the 
real 3D-width of the shower varies significantly with the zenith angle, a new parameter called 
the reduced 3D-width \6\ was introduced, which is the ratio of the 3D-width to the slant depth at 
shower maximum, both being expressed in g/cm^ at the altitude of shower maximum. This last 
variable was found to be almost independent of the zenith angle. 

The other crucial observable quantity examined, the logarithm of the number of Cherenkov 
photons, is related directly to the energy measurement. As this is the single parameter in the 
3D-model that has to be calibrated from simulations, its verification is essential. Fig. [12] (left 
panel) and Table |6] show a good agreement between data and simulations on the distribution of 
the logarithm of the number A^^. Cherenkov photons, when the optical efficiency derived from 
muon rmgs is used in both cases. As a matter of fact, using the same optical efficiency 
in simulation and in data analysis is not sufficient to guarantee a good agreement between the 
distributions of In(A^c)- This is shown in Fig. [12] (right panel), in which this common efficiency 
has been increased by 10% with respect to that given by muon rings resulting in a manifest shift 
of the ln(A^c) distributions. Therefore, the agreement shown in Fig. [12] (left panel) and Table [6] 
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Figure 12: Comparison of the number of reconstructed Cherenkov photons for the data (red points with error bars) 
and 3D-Model simulations (black histograms). Left panel: the optical efficiencies have been measured from muon ring 
images. Right panel: the optical efficiency has been increased by 10%. 



implies that the optical efficiency used in the simulations (both in the event generation and in its 
reconstruction) matches the real one. In contrast to the previously studied variables, the energy 
measurement seems to be the one most affected by variations in the experimental conditions. In 
particular, its sensitivity to the relative optical efficiency is critical and emphasises the need for a 
careful calibration. 



Observable X 


100 X ^ 


A(HWHM) 
Bin width 


Number of pixels per telescope 


-2.82 + 3.00 


0.62 


Charge in the 3"* brightest pixel 


2.60+ 1.23 


0.17 


Altitude 


-0.22 ± 0.37 


-0.51 


Slant depth 


-1.11+0.55 


-0.06 


3D-width 


-1.53 + 0.72 


0.31 


Reduced 3D-width 


-0.91+0.67 


0.33 




0.18 + 0.10 


-0.47 



Table 6: Comparison between simulated and real distributions shown in Fig. \9\ \10\n7\ and \12\ (left panel). For each 
observable X, whose mean values are Xj (real data) andXs (simulations), the relative bias 2(Xj -X,)/{Xj + Xs) = AX/X 
is given in percent units in the second column. Each distribution has also been characterised by its half-width half- 
maximum (HWHMj for data, HWHM, for simulations). In the last columnn, the difference A{HWHM) = HWRMj - 
HWHM, is expressed as a fraction of the bin width. 



3.4. Quality criterion, gamma-ray selection efficiency and background rejection 

Finally, the performance of the method determined from simulations and discussed in sec- 
tion l2.3l can be checked on the basis of real data from the flare of PKS 2155-304. The experimen- 
tal distribution of the variable u - ln(23D/2reai) defined for triggered telescopes in section IZTl 
shown by points with error bars in Fig. [13] (left panel) for an average zenith angle of 18° and 
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l"(QiD/Qrea,) ln(Q3D) 



Figure 13: Left panel: comparison of the distribution of the variable u = \n{QT,o I QmA) for the data (points with error 
bars) for an average zenith angle of 18° and nj = 1, to the curve deduced from formula JT} and Table\l\ Right panel: 
comparison of the distribution of the variable IniQio) for the data (points with error bars) in the same conditions, to the 
Gaussian distribution whose parameters are deduced from Table\2\ 



nj - 2, shows a good agreement with the curve deduced from formula ([TJ and Table [T] In the 
right panel of Fig. [13] the experimental distribution of \vi{Q^D) for non triggered telescopes is 
found to be in good agreement with the Gaussian whose parameters are deduced from Table |2l 
It should be noted that both curves were obtained from simulations based on an energy spectrum 
with photon index 2.2, whereas the corresponding index for PKS 2155-304 during the flaring 
period is 3.3 ifisll . This results confirms, on the basis of real data, that the distributions of the 
new variables u and v are not sensitive to the spectral index of the gamma-ray source under study, 
as already checked from simulations. The hadronic rejection factors at 18° zenith angle obtained 
from PKS 2155-304 data (5085 + 343 for the 3D standard analysis and 7935 + 668 for the new 
one), are found in good agreement with the figures given in Table |4] obtained from simulations 
for point-like sources. Only the ratio of gamma-ray selection efficiencies in the new and in the 
standard 3D analyses respectively can be checked experimentally. At 18° zenith angle, this ratio, 
found to be 88% in PKS 2155-304 data, is very close to the value of 87% obtained from simu- 
lations and derived from Table |4] The preceding agreement between data and simulations, both 
on hadronic rejection factors and on gamma-ray selection efficiencies, validates the estimation 
of the gain in sensitivity achieved by the new 3D analysis given by Tables[3]and|4] 

4. Conclusion 

The reconstruction method based on the 3D model of electromagnetic showers has been im- 
proved by an additional consistency check: namely, the modelled Cherenkov photosphere must 
reasonably reproduce the global light yield of all triggered telescopes and be compatible with a 
low amount of light for non-triggered ones. The new criteria involved in this consistency check 
are quasi-insensitive to the energy spectrum of gamma-rays used in the simulations. This was 
already the case for the gamma-ray/hadron discrimination criteria described in the 3D standard 
analysis ^ and essentially based on well-known properties of electromagnetic showers: rota- 
tional symmetry and lateral spread of the Cherenkov photosphere at shower maximum. This 
makes this method quite different from those based on machine learning algorithms or training 
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procedures (e.g. whose sensitivity is optimised for a given gamma-ray spectrum. The gain 
in sensitivity with respect to the 3D standard analysis corresponds to that of an observation time 
increased by 50% for an extended source at low zenith angles, the hadronic background being 
lowered by a factor of 2. This contamination is particularly reduced for the sample of events 
triggering only two telescopes, i.e. for the majority of events. Finally, the almost pure gamma- 
ray beam provided by the exceptional flare of the blazar PKS 2155-304 in July 2006 allowed 
to check the simulation of the H.E.S.S. experiment at the level of a few percent and to validate 
the gamma-ray energy measurement based on the 3D analysis. The present method, which takes 
account of the correlations between images observed by different telescopes and allows a rather 
fast processing of the events, will be particularly well suited to the exploitation of forthcoming 
large Cherenkov arrays. 



5. Acknowledgments 

We thank Prof. W. Hofmann, spokesman of the H.E.S.S. Collaboration and Prof. G. Fontaine, 
chairman of the Collaboration board, for allowing us to use H.E.S.S. data in this publication. We 
are grateful to Dr. M. de Naurois for carefully reading the manuscript and for providing us with 
very useful suggestions. Finally, our thanks go to all the members of the H.E.S.S. Collaboration 
for their technical support and for many stimulating discussions. 



References 

[1] M. Punch, in "Towards a Network of Atmospheric Cherenkov Detectors VII", Palaiseau, France (2005), p. 379, 

B. Degrange et G. Fontaine ed. 
[2] M. Mori, in "Towai'ds a Network of Atmospheric Cherenkov Detectors VII", Palaiseau, France (2005), p. 19, 

B. Degrange et G. Fontaine ed. 
[3] K. Bemlohr, et al., Astropart. Phys. 20 (2003) 111. 

[4] T.C. Weekes, in "Towards a Network of Atmospheric Cherenkov Detectors VII", Palaiseau, France (2005), p. 3, 
B. Degrange et G. Fontaine ed. 

[5] M. Teshima, in "Towards a Network of Atmospheric Cherenkov Detectors VII", Palaiseau, France (2005), p. 373, 
B. Degrange et G. Fontaine ed. 

[6] M. Lemoine-Goumard, B. Degrange and M. Tluczykont, Astropart. Phys., 25 (2006) 195. 

[7] S. Funk, et al., Astropart. Phys. 22 (2004) 285. 

[8] D. Berge, S. Funk and J. Hinton, A&A 466 (2007) 1219. 

[9] F Aharonian, et al., A&A 455 (2006) 461 . 
[10] F Aharonian, et al., ApJ 661 (2007) 236. 
[II] G. Vacand, et al., A.stropart. Phys. 2 (1994) I. 
[12] F. Aharonian, et al., Astropart. Phys. 22 (2004) 109. 
[13] F Aharonian, et al., ApJ 664 (2007) L7I. 

[14] J. Guy, these de doctorat, Universite Paris VI (2003); this program is a modified version of that described in 

M. P Kertzman and G. H. Sembroski, Nucl. Instr Meth. Phys. Res. A, 343 (1994) 629. 
[15] J. Guy, these de doctorat, Universite Paris VI (2003). 
[16] F Aharonian, et al., A&A 457 (2006) 899. 

[17] T.K. Gaisser, in Cosmic Rays and Particle Physics, Cambridge University Press (1990), p. 245. 

[18] M. de Naurois, in "Towards a Network of Atmospheric Cherenkov Detectors VII", Palaiseau, France (2005), p. 

149, B. Degrange et G. Fontaine ed. 
[19] S. Ohm, C. van Eldik and K. Egberts, Astropart. Phys. submitted. 

F. Dubois, G. Lamanna and A. Jacholkowska, Astropart. Phys. submitted. 



17 



